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Abstract 

H , 

^ Cy ^ We investigate the dynamics of a flat isotropic brane Universe with 

two-component matter source: perfect fluid with the equation of state 
p = (7 — l)p and a scalar field with a power-law potential V ^ 4)°". The 
index a can be either positive or negative. We describe solutions for which 
the scalar field energy density scales as a power-law of the scale factor (so 
called scaling solutions). In the nonstandard brane regime when the brane 
is driven by energy density square term these solutions are rather different 
from their analogs in the standard cosmology. A particular attention 
is paid to the inverse square potential. Its dynamical properties in the 
. nonstandard brane regime are in some sense analogous to those of the 

exponential potential in the standard cosmology. Stability analysis of 
fSJ \ the scaling solutions are provided. We also describe solutions existing in 

. regions of the parameter space where the scaling solutions are unstable or 

' do not exist. 

7'. 

1 Introduction 

Cosmological evolution of a brane world has become a matter of intense in- 
Jv>( \ vestigation after the famous paper of L. RandaU and R. Sundrum 1 . Due to 

}_( ' additional terms, which appear in effective Einstein equations on a brane, the 

brane cosmology has some important differences from the standard cosmological 
scenario. 

The field equations on the brane are |21IS| 

C^j/ — K Tfiii + h S^i, fpi/. (1) 

^ 3/2 

Here k — is the fundamental 5-dimensional gravitational constant, k = 

m^^ is the effective 4-dimensional gravitation constant on the brane. Bulk 
corrections to the Einstein equations on the brane are of two forms: there are 
quadratic energy-momentum corrections via the tensor S^i, and nonlocal effects 
from the free gravitational field in the bulk 5^^. The nonlocal corrections can 
be decomposed into a scalar, vector and tensor parts 0- On FRW brane there 
is only a scalar part U which decays during the expansion of the Universe as 
U ~ a^^ where a is the scale factor. This means that this term behaves exactly 
as a radiation fiuid, and for this reason it is often called as "dark radiation" 
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despite its pure geometrical nature. It should be noted, however, that "dark 
radiation" can be either positive or negative, and in the latter case it gives 
rise to some cosmological solutions, impossible in the standard cosmology (for 
example, the collapse of a flat brane Universe U). 

More important differences between standard and brane cosmologies appear 
at early stage of the Universe evolution, where quadratic corrections propor- 
tional to S^j^u dominate over the standard terms proportional to T^^. The rate 
of expansion of the Universe changes with respect to matter density, this leads 
to interesting peculiarities of the brane scenario. 

The case of matter source in the form of perfect fluid with the equation of 
state p = (7 — l)p have been extensively investigated during last three years, 
see, for example, 

The next step is to consider a scalar field on a brane. Even in the standard 
cosmology the dynamics of Universe filled by a scalar field is much more com- 
phcated in comparison to the perfect fluid case [Jj. The efi^ective 7 of the fleld 
can vary in the whole possible interval 7 S (0,2). This property, if the scalar 
field potential is not too steep, leads to existence of inflation and, in the case 
of positive spatial curvature, to the possibility of bouncing solutions and dy- 
namical chaos. However, if the potential is steep enough, the situation changes 
and this transition occurs for potentials as steep as exponential: l^(^) ~ exp 
with large enough A can not provide inflation [5] and bounces It already 
have been noticed that in brane cosmology we have a similar picture but with 
inverse square potential V{<t)) = B<j>~^ playing the role of exponential potentials 
in the standard cosmology [IT)] . 

In the present paper we consider a brane Universe filled by a scalar fleld 
and a perfect fluid. In the ordinary cosmology exponential potentials in this 
situation have again a particular property - in this case a tracker solution (the 
energy density of a scalar field tracks exactly the energy density of an ordinary 
matter) exists ^] ^] . Not surprising that a tracker solution on a brane have 
been found for the inverse square potential. It was done by K.-I. Maeda ^3] in 
the particular case of ordinary matter in the form of radiation fluid (7 — 4/3). 
We present a general form of brane tracker solution for an arbitrary 7. We 
also study other power-law potentials which give us brane scaling solutions - 
solutions with scalar field energy density scales as some power of a scale factor, 
when the perfect fluid, being the dominant component, has an energy density 
scaling as a possibly different power (we use the terminology of 52)- Such 
potentials are often used in modern quintessence models ^3 1141 1151 III)] to 
explain observable acceleration of our Universe. 

Our analysis can be useful for understanding evolution of a quintessence field 
before the radiation-dominated stage, for example, during inflation. Though for 
the full treatment of this problem it is necessary to consider a more complicated 
case of a two-field dynamics (quintessence and inflaton fields), we can use the 
approximation of a perfect fluid if 7 of the matter varies slowly enough. 

In the present paper we consider only a flat isotropic brane. Shear and cur- 
vature terms can modify the scaling solutions, in the standard cosmology it was 
shown in A two-component (a scalar fleld and a perfect fluid) brane 
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dynamics with nonzero shear and curvature has aheady become a matter of in- 
vestigations |19[I20| . We leave the problem of scaling solutions on an anisotropic 
brane to a future work. 



2 Tracker solutions 
2.1 Basic equations 

We start with system of equations describing a scalar field on a flat isotropic 
brane There are the Raychaudhuri equation 



H = -H'_^£(, + 3rt-^(2.H-W-|i, (2) 

the Friedmann equation 



the Klein-Gordon equation for the scalar field 

^ + 3H:j> + V'{(t)) ^0, (4) 
and the evolution equation for the " dark radiation" 

U + \hU = 0. (5) 

Here A is the brane tension. The overall pressure and energy density in (2) - 
(3) are the sums 

P=Pm+P4>, P = Pm+P4,, 

where the perfect fluid equation of state Pm — (cim ^ ^)Pm, 7m = const and the 
scalar field equation of state p^f, — (7^ — l)p0 with 7^ changing with time. The 
pressure and the energy density of the scalar field are given by 



Using Eq. (4) we can write the dynamical equation for 7^ as 

V 

70 = 3iJ70(70 - 2) - 270^. (6) 

In this section we study the particular potential V{(j)) — B(f)^^ , for this case the 
Eq.(6) becomes 

70 = 3i?70(70-2) + ^^P0(2-70)'/^ (7) 
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We would like to notice an interesting difference between scalar field dynam- 
ics in ordinary cosmology and on a brane. In standard cosmology a constant 
coefficient in field potential can be absorbed by rescaling scale factor and time 
without changing a general dynamical picture. In brane cosmology it is not the 
case, and coefficient B plays an important role, as we will see later. 

After introducing the deceleration parameter q as 



dimensionless variables 



and a new time variable r such that ^ = we can rewrite our system in the 
form 

q = ^^^pm(37„i - 2) + ir2p0(370 - 2) + ^\r,i{ilm - 1) + 



f^A0(370 - 1) + V^^A^n^^(3(7m + 70) " 2), (8) 



1 = Upm + + U + r^Am + ^^A0 + 2 -\/TT\mT^A0 , (9) 

n'p„,^npra['2{l + q)^ilral ^'p^ ^ np^[2{\ + q) - i^^], (10) 

f^Am -f^Am[2(I + (7)-67m], f^^^ - A0 [2(1 + ?) - 67^] , (II) 

U' = 2U{q~l), (12) 



7; = 3747. - 2) + ^l^^^ml{2 - ^,f\ (13) 

Here the prime denoted derivative with respect to r. 

Note, that in the case of several matter sources nonlinear crossing terms 
appear in the r.h.s. of Eqs. (8) - (9). 

In this paper we study only a pure nonstandard regime when the contribution 
from linear terms proportional to Tp^ in (1) is negligible in comparison to the 
quadratic energy- momentum terms proportional to Sp^. Therefore we put ^Ipm 
and rip0 to zero. For the potential V{4)) = B(j>^^ the Hubble parameter H 
decouples (in standard regime this property belongs to an exponential potential 
[22j) and we have a compact phase space unless a "dark energy" density U < 0. 
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2.2 Equilibrium points 

Now we list equilibrium points for the system (8) - (13) which can be stable in 
future direction with corresponding results of their stability analysis. To sim- 
plify the formulae we use the dimensionless coefficient b = B/m^. Remember, 
that and ilp^ are set to zero. As for the full system (8) - (13) the equilib- 
rium points obtained correspond to an intermediate brane behavior, when some 
asymptotic regime have already been established, but the brane is still driven 
by the quadratic terms in (1). First, consider the case of ?7 = 0. 



Stable for 7„ > j:^- 

II 

70 = Im, "A0 = 7, "Am — 



8(2 -7m)' 8(2 -7„) ^2(2 -7„)' 

Exists and stable for 7„i < 

The point I is the regime of scalar field dominance. We see that depending 
on b the scalar field density can mimic an ordinary matter in the whole interval 
of possible values 7^ £ (0,2). Naturally, this point is stable for 7^ > 7^. The 
point II represents the brane tracker solution with flx^/flx^ — const, which 
was obtained in for the particular case of 7^ = 4/3. 

The case of [/ > leads to the following changes. The point I disappears 
for b > 16, the point II becomes unstable for 7„i > 2/3. In these regions a new 
stable point appears: 

III 

70 = 2/3, ^A0 = Y^, ^ = 1-1^' ^^Am = 0. 

Exists for b < 16, stable for 7^ > 2/3. 

This is an analog of tracker solution when the " dark energy" plays the role of 
matter. The ratio fl\cj,/U remains constant. Of course, for 2/3 < 7m < 4/3, 
this situation (" dark energy" dominates the ordinary matter) changes after the 
low-energy stage begins and flpm becomes important. Zones of stability for the 
future direction in the plane (6, 7^) for the case of [/ > are plotted in Fig.l. 

If U is negative, the phase space becomes incompact. Physically this means 
that a recoUapse of the brane Universe is possible. Our numerical studies in- 
dicate that in zones of stability of points I and II the recollapse occurs for 
large enough initial \U\, while for small \U\ the Universe tends asymptotically 
to corresponding stable points I or II with U — 0. On the contrary, in the 
region III there are no stable equilibrium points, corresponding to initial neg- 
ative U. We have confirmed numerically, that all trajectories from the region 
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Figure 1: Stability zones of the equilibrium points for the case U > on the 
plane (&, 7m) where b = B/m^. 

b < 16,7m > 2/3, U < recollaps. It is worth to remind, however, that all this 
results are obtained for a pure nonstandard regime. When the full system of 
equations (8) - (13) is considered, some trajectories may still not recollaps until 
transition to the standard regime occurs, and, if 2/3 < 7m < 4/3 the recollapse 
will never occur. 

3 General power-law potentials 
3.1 Scaling solutions 

In this section we describe solutions with scalar field energy density scales ex- 
actly as some power (in general, different from the fluids one) of the scale factor. 
The case of scalar field dominance have been described in , therefore we con- 
centrate our effort to the fluid dominance case. Our work generalize PB| to fluids 
with 7 ^ 4/3. 

It is worth to note that these solutions are not future attractor in the sense 
of the previous section even if we do not consider the transition to the standard 
regime. Indeed, the dynamical evolution may break down the fluid dominance. 
Nevertheless, it is an interesting example of transient behavior. This behav- 
ior is followed by a regime of a scalar fleld dominance. Late-time scalar field 
dominance in the standard regime leads to inflation 12 and provides a possible 
explanation of accelerated expansion of our present Universe. 

We will follow the method of where a similar class of solutions have 
been found in the standard cosmology. Suppose that the scalar field energy 
density behaves as ^ where a is the scale factor. Since the Klein- 

Gordon equation is not modified on a brane, this property, as in the standard 
case, requires that the ratio of a scalar field kinetic energy density and a total 
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scalar field energy density remains constant |12| : 



^ = I (14) 

The perfect fluid with equation of state p = (7 — l)p has the energy density 
p ~ a~™ where m — 87. In the pure nonstandard brane regime the fluid 
dominance leads to 

a-ii/". (15) 

All our considerations remain valid also for the case of the " dark energy" dom- 
inance (in this case, of course, > 0), because the "dark energy" behaves as a 
"normal" matter with m = 4. Due to (15), Eq.(4) becomes 

I 6 l-dV 

= T'^TJT- (16) 

m t aq) 

Eq. (14) gives us 

Integrating this equation we obtain 

(j) = Af-"/^"". (17) 

Substituting (17) into (16) and solving for V{(j)) we find the potentials which 
allow the scaling behavior in the nonstandard brane regime 

n^)^A^- ^(^("-^)"rV ", (18) 

am[a — 2)^ 

where 

2n 



o • (19) 

n~-2m ^ ^ 

It means that, as in the standard cosmology, scaling solutions exist for power- 
law potentials V{(l)) ^ cf)"-. If a is given, the scalar field energy density is 
proportional to a~", where 

n = m. (20) 

a — 2 

Though this expression differs only for the factor 2 in the numerator from its 
analog in the standard cosmology the properties of brane scaling solutions 
are rather different. Let us begin with the case of negative a. In the standard 
scenario the energy density of a scalar field with a < always falls slower than 
the density of a perfect fluid while in the brane cosmology it falls slower for 
a > —2 and faster for a < —2. Second, the positivity of the potential in (18) 
requires that a(3 — m) < 6, which indicates that for m > 3 (a matter with a 
positive pressure) and 

« < (21) 
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the scaling solution does not exist. On the contrary, in the standard cosmology 
the scaling solution exists for an arbitrary negative a independently on m. 

Consider now the case of positive a. The positivity of (18) leads to a(3 — 
m) > 6. It means that for m > 3 scaling solutions do not exist for any a. In 
particular, they do not exist for the case of 7 = 4/3, investigated in [531 ■ 
the other side, for to < 3 the situation is similar to the standard cosmology - 
scaling solution exists for sufficiently steep potential with a satisfying in the 
brane case 

6 , ^ 

" > (22) 



3 — TO 



3.2 Stability analysis 



We have to analyze the stability of the scaling solutions (17) with a power-law 
potential (18). To do this we use new variables 

T = e*, u{t) = ^ , (23) 

where ipo{t) is the exact solution given by (17) Then the equation for scalar 
field becomes 

\2 — a m J a — 2 \a — 2 m J ^ ' 

where the prime denotes the derivative with respect to r. This second-order 
differential equation can be divided into the autonomous system 

Scaling solution corresponds to a critical point (u,p) — (1,0). Linearizing 
51 about this point wc find the eigenvalues of these coupled equations 




3 \^ 2a 6 



^^■^ 2(a-2) 2TO^\/v2(a-2) 2m j a - 2 to' 

Then the condition for stability is the negativity of the real part of both eigen- 
values. 

For a < after analysis wc find that scaling solution is a stable attractor 
for all region of existence (21). And for a > there is an additional inequality 

2™ + 6 

a > , 27 

3 — TO 

which excludes a narrow domain from (22). For example, for ~ 0^ the 
scaling solution is stable if m < 1, or, equivalently, 7^ < 1/3, this is exactly 
the condition for a brane inflation driven by an ordinary matter. For 1/ ^ 
we have to < 3/2 (or 7^ < 1/2). 

Zones of stability of the scaling solutions on the plane (a, to) are plotted in 
Fig.2. 
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Figure 2: Stability zones of the scaling solutions (dark). In the white zone with 
negative a the kinetic-term dominated solution is stable, in white zone with 
positive a the scalar field infinitely oscillates. 

3.3 Other solutions 

In the end of this section we briefly describe solutions, existing in the regions of 
the plane (a,m) where scaling solutions do not exist. 

For positive a with even n we have infinite damping oscillations of the scalar 
field, similar to the standard case. Some asymptotic formulae for potentials 
V = ^m?(fp and ^A0'* have been given in |23j . 

For negative a the other type of solution exists. In this kinetic-term dom- 
inated solution the contribution of the potential for the energy density of the 
scalar field is negligible in comparison to the contribution of the kinetic term. 
Its explicit form is 



(28) 



where (j)Q and io are integration constants. This solution does not exist for 
m < 3. In this case the scaling solution is stable for all negative a. For m = 4 
we recover the known behavior '-^ t^/"* ^23|. 

4 Conclusions 



We have investigated a scalar field dynamics with a power-law potential on a 
brane filled by a perfect fluid with equation of state Pm = {jm — ^)pm in a regime 
when quadratic energy-momentum correction dominates. The tracker solution 
when the ratio of scalar field and fluid energy densities remains constant is found 
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for the scalar field potential V{(j)) = B(f)~'^ . Two other future asymptotic for this 
case are described and zones of their stability in the plane (i?,7m) are found. 

We also describe scalar field evolution with V{(j)) ~ in a brane Universe, 
dominated by perfect fluid energy density. We have found scaling solution, 
analogous to known in the standard cosmology. The main differences between 
standard and brane cosmology are: 

• For negative a in the standard cosmology the scaling solutions exist for 
arbitrary and a, whereas in the brane cosmology they do not exist for 
7m > 1 and large enough |a|, which do not satisfy (21). 

• For positive a in the standard cosmology the scaling solutions exist for 
any 7^ if a exceeds some value, depending on 7„j. In the brane cosmology 
the scaling solutions exist only for 7„i < 1 and large enough a, satisfying 
(27). They are absent for 7,„ > 1. 

It is also interesting, that in the case of negative a the scalar field energy 
density may fall more rapidly or more slowly than the fluid density, depending 
on the sign of a + 2, whereas in the standard scaling solution with negative a 
the ratio of scalar field and fluid energy densities always increases. This was 
discovered in the particular case of radiation fluid in ^21 1 where it was argued 
that this fact may be important for quintessence scenario. In our paper we 
confirm that this property do not depend on the state equation of the perfect 
fluid. 
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